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» I ' The random-phase-approximation semiclassical scheme for description of plasmon excitations in 

, large metallic nanospheres (with radius 10-100 nm) is developed for a case of presence of dynamical 

electric field. The spectrum of plasmons in metallic nanosphere is determined including both surface 
and volume type excitations and their mutual connections. It is demonstrated that only surface 
QQ ] plasmons of dipole type can be excited by a homogeneous dynamical electric field. The Lorentz 

^SJ ■ friction due to irradiation of e-m energy by plasmon oscillations is analysed with respect to the 

sphere dimension. The resulting shift of resonance frequency due to plasmon damping is compared 
with experimental data for various sphere radii. Collective of wave-type oscillations of surface 
plasmons in long chains of metallic spheres are described. The undamped region of propagation of 
plasmon waves along the chain is found in agreement with some previous numerical simulations. 
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I. INTRODUCTION 



Experimental and theoretical invesitigations of plasmon excitations in metallic nanocrystals rapidly grew mainly 
due to perspectives of possible applications in photovoltaics and microelectronics. A significant enhancement of 
' absorption of the incident light in photodiode-systems with active surface covered with nano-dimension metallic 
' particles (of Au, Ag or Cu) with planar density ^ 10^/cm^ was observed^"— . This is due to a mediating role in light 
(—| ■ energy transport played by surface plasmon oscillations in metallic nano-compounds on semiconductor surface. These 
Q findings are of practical importance towards enhancement of solar cell efficiency especially for thin film cell technology. 
O . Hybridized states of the surface plasmons and photons result in plasmon-polaritons^ which are of high importance 
for applications in photonics and microelectronics^^, in particular for transportation of energy in metallic modified 
^-H • structures in nano-scal o^^'^^ . 

^ ' Surface plasmons in nanoparticles were widely investigated since their classical description by Mie^^. Many partic- 
ular studies, including numerical modelling of multi-electron clusters, have been carried ou^^—. They were develop- 
ments of Kohn-Sham attitude in form of LDA (Local Density Approach) or TDLDA (Time Dependent LDA)^ii^iii'^ 
addressed, however, to small metallic clusters, up to ca. 200 electrons (as limetted by severe numerical constraints). 
The random phase approximation (RPA) was formulated^^ for description of volume plasmons in bulk metals and 
utilised also for confined geometry mainly in numerical or semi- numerical manner^^. Usually, in these analyses the 
jellium model was assumed for description of positive ion background in metal and dynamics was addressed to elec- 
tron system onlyi^"— , and such an attitude is preferable for clusters of simple metals, including noble metals (also 
transition and alkali ones). 

In the present paper we generalise the bulk RPA descriptioniS, using semiclassical approach, for a large metallic 

nanosphere (with radius of several tens nm, and with 10^ 10^ electrons) in an all analytical calculus version—. 

The plasmon oscillations of compresional and traslational type, resulting in excitations inside the sphere and on its 
\ surface, respectively, are analysed and referred to volume and surface plasmons. Damping effects of plasmons via 
electron scattering processes and radiation losses are included, the latter ones, via Lorentz friction force. The shift 
of the resonance frequency of dipole-type surface plasmons (only plasmons induced by homogeneous time-dependent 
electric field), due to damping phenomena, is compared with the experimental data for various nanosphere radii. 
Collective surface dipole-type plasmon oscillations in the linear chain of metallic nanospheres are analysed and wave- 
type plasmon modes are described. A coupling in near field regime between oscillating dipoles of surface plasmons 
together with retardation effects of energy irradiation lead to a possibility of undamped propagation of plasmon waves 
along the chain in the experimentally realistic region of parameters (of separation of spheres in the chain and their 
radii). These effect would be of particular significance for by plasmon arranged transport of energy along metallic 
chains for application in nanoclcctronics. 

The paper is organised as follows. In the next section the standard RPA theory in quasiclassical limit, is generalised 
for the confined system of spherical shape. The resulting equations for volume and surface plasmons are solved in 
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the following section (with particularities of calculus shifted to the Appendix). The next section contains description 
of the Lorentz friction for surface plasmons oscillations of the dipole-type. The analysis of the collective wave-type 
surface plasmon oscillations in the chain of metallic nanospheres is presented in the last section. Besides the theoretical 
model the comparison of the characteristic nano-scale plasmon behaviour with available experimental data, including 
own measurements, is presented. 



II. RPA APPROACH TO ELECTRON EXCITATIONS IN METALLIC NANOSPHERE 
A. Derivation of RPA equation for local electron density in a confined spherical geometry 



Let us consider a metallic sphere with a radius a located in the vacuum, e = 1, /i = 1 and in the presence of 
dynamical electric field (magnetic field is assumed to be zero). We will consider collective electrons in the metallic 
material. The model jelliurr}^^^^ is assumed in order to account for screening background of positive ions in the form 
of static uniformly distributed over the sphere positive charge: 



ne(r) = neQ{a - r). 



(1) 



where Ue = N^/V and ne\e\ is the averaged positive charge density, N,, the number of collective electrons in the sphere, 

V = '^^^ the sphere volume, and Q is the Heaviside step-function. Neglecting the ion dynamics within jellium model, 
which is adopted in particular for description of simple metals, e.g. noble, transition and alkali metals, we deal with 
the Hamiltonian for collective electrons, 
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ne{ro)d^ro 
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eip{R + rj,t) 
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(2) 



where Vj and m are the position (with respect to the dot center) and the mass of the j electron, R is the position of 
the metallic sphere center, AE represents electrostatic energy contribution from the ion 'jellium', ip{r,t) is the scalar 
potential of the external electric field. The corresponding electric field E{R + rj,t) = —gradj(p{R + rj,t). Assuming 
that space-dependent variation of E is weak on the scale of the sphere radius a then E{R + rj,t) ~ E{R, t), i.e the 
electric field is homogeneous over the sphere (it holds for \R\ ^ a). Then (p{R + rj,t) ~ —E{R, t) ■ R+Lpi{R + rj,t), 
where (pi{R + rj,t) ~ —Vj ■ E{R, t). Hence, one can rewrite the Hamiltonian ([2]) in the form: 



He = H'^- eNE{R,t) ■ R, 



where 
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the corresponding wave function can be represented as. 



+ eip{rj,t) 
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AE, 



<\,t j-\ n,i I j-\ f E(R,t)-Rdt 



with ih^ = Him 



re = (r-i,r2, ...,rAr). 
A local electron density can be written as followgl^.: 



p(r, t) =< vl/(re, t)\Y, - r,)\mire,t) >=< m'ir,,t)\ ^ 6{r - r,)|*'(r„ t) >, 

j j 

with the Fourier picture: 
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where the 'operator' p{k) = 



3 



Using the above notation one can rewrite H'^ in the following form, in analogy to the bulk casi 
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+ ^ J d^k^,{k,t) (p+{k) + m) + ^ I d'k^ p+{k)p{k)^N, 
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where: ne{k) — J (Prne{r)e 
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Utilizing this form of the electron Hamiltonian one can write the secod time-derivative of p{k): 

d^p{k,t) _ 1 



dt^ 



(9) 



which resolves itself into the equation: 



d^Spjk.t) _ 



l^Jd^qfie{k-q)!^Sp{q) - J d^qne{k - q){k ■ q)<fi{q,t) 



m2 

J d\5p[k -q){k- (q, t)~^! d^qSpik - q)^Sp{q), 



(10) 



where Sp{k) — p{k) — ne{k) is the 'operator' of local electron density fluctuations beyond the uniform distribution. 
Taking into account that: Sp{k,t) =< "ii' {t)\Sp{k)\^' (t) >= p{k,t) - fieik) we find: 



d^Sp(k,t) _ 



> 



J d^qfieik - q)^Sp{q, t) - ;J;;3 / d^qne{k - q){k ■ q)(pi{q, t) 
1^ J d\5~p{k ~ q, t){k ■ t) - ^ J d\^ < ^'\5p{k ~ q)Sp{qW >, 



(11) 



One can simplify the above equation upon the assumption that Sp{r, t) = J e^^''^5p{k^ t)d^k only weakly varies on 
the interatomic scale, and hence three components of the first term in right-hand-side of Eq. can be estimated as: 

k'^VpSp{k), k'^vp /kT5p{k) and k^Vp/k'^Si5{k), respectively, with l/kx the Thomas-Fermi radius^^, fcy = /eTrn^e^ 



jpu/jyruj, n, u ^ / ri^ u fjyru j cxiiu n, u p / rirpi, ^ ^. . ^ , ^ ^, , ..^ , ..^ 

ep the Fermi energy, and vp the Fermi velocity. Thus the contribution of the second and the third components 
of the first term can be neglected in comparison to the first component. Small and thus negligible is also the last 
term in right- hand- side of Eq. (Ilip . as it involves a product of two dp (which we assumed small Sp/ug « 1). 
This approach corresponds to random-phase-approximation (RPA) attitude formulated for bulk metali^i^i (note that 
(5/5(0) = and the coherent RPA contribution of interaction is comprised by the second term in Eq. The last 

but one term in Eq. (1111) can also be omitted if one confines it to linear terms with respect to dp and ipi. Next, due 

to spherical symmetry, < ^ e^*'" ''^^ (fc • Vj)^ l^*' >~ |^ < ^''| ^^2^1^' >■ Performing the inverse 



Fourier transform, Eq. (|lip attains finally the form: 



d,t 



-gv {e(a ~ r)V / d\,^^5p{r,,t)] + {e(a - r)V^i (n, i)} . 



(12) 



According to Thomas-Fermi approximation the RPA averaged kinetic energy can be represented as follows: 
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(13) 



Taking then into account the above approximation and that V0(a — r) = —^5{a — r) = lim£_j.o 5{a + e — r) as 
well as that (pi{R, r, t) — —r ■ E{R, t), one can rewrite Eq. ([T^ in the following manner: 

^^^[l^^'Spir,t)-u:lSpir,t)]e{a-r) 

V { [lepn, + ep6p{r, t)] r ^(a + e - r) } - [| ^ ^ V<5p(r , t) + ^^\/J d^n p^'5p(ri ,t) + ^^-E{R, t)] S{a + e-r). 
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In the above formula ujp is the bulk plasmon frequency, Wp = '^^^'^ , and S{a + e — r) = linij_j.o (5(a + e — r). The 
solution of Eq. can be decomposed into two parts related to the domain: 

Sp{r,t)=llPf^'} (15) 
\ 5p2{'r,t), for r > a, (r ^ a+), ^ ' 

corresponding to the volume and surface excitations, respectively. These two parts of local electron density fluctuations 
satisfy the equations: 

^!^ = ^^V^^p,(r.,t)-.>,(r,.), (16) 

and 

- -alrV { [f e^n. + .,8p,[r^ t)\ r<5(a + . - r)} 
- |^^V<5p2(r,t) + It^V/d^n^ (<5pi(ri,t)e(a - n) + 5p2(ri, t)e(ri - a)) + • £;(i2,i)] <5(a + e- r). 

(17) 

It is clear from Eq. (|16p that the volume plasmons are independent of surface plasmons. However, surface plasmons 
can be excited by volume plasmons due to the last term in Eq. (jl7p , which expresses a coupling between surface and 
volume plasmons in the metallic nanosphere within RPA semiclassical picture. It is in fact a surface tail of volume 
compressional-type excitations, while surface traslational-type exctations have no a volume tail. 

In a dielectric medium in which the metallic sphere can be embedded, the electrons on the surface interact with 
forces e (dielectric susceptibility constant) times weaker in comparison to electrons inside the sphere. To account for 
it, one substitutes Eqs (|T6| and (flT)) with the following ones: 

dt^ ' - V^^pi(r,t)-^>i(r,^), (18) 

and 

= { [f^Fne + epSp,ir, t)] ^S{a + e - r)} 

- I^^VSp2{r,t) + grv/d^ri^ (Vi(n,t)e(a - n) + i<5p2(ri,i)e(ri - a)) + ■ EiR,t)\ S{a + e- r). 

(19) 

Let us also assume that both volume and surface plasmon oscillations are damped with the time ratio tq which can 
be phenomenologically accounted for via the additional term, — :^^^^f^) to the right- hand- side of above equations. 



They attain the form: 

d^6pi{r,t) 2 dSp{r,t) 2 ep 



V^Spiir,t)-u;lSp,ir,t), (20) 



dt^ tq dt 3 TO 

and 

a^S^ + ^-i,^ { [hFU. + epSp^ir, t)] r^5{a + e - r)] 

- |f^^V<Sp2(r,t) -I- S^V/d^ri^ (5pi(ri,t)e(a - n) + \5p^{r^,t)Q{n - a)) + . E{R,t)\ 5{a + e- r). 

(21) 

From Eqs ([^0]) and (l^lT) it is noticeable that the homogeneous electric field does not excite the volume-type plasmon 
oscillations but only contributes to surface plasmons. 



B. Solution of RPA equations: volume and surface plasmons frequencies 

Eqs (|20|) and (|2T|) can be solved upon imposing the boundary and symmetry conditions — cf. Appendix [A] Let us 
write the both parts of the electron fluctuation in the following manner: 

5 Pi ir,t)= rie [/i (r) + F{r, t)] , for r < a, ,^2) 
Sp2ir,t) =nef2ir) +a{n,t)S{r + a), for r > a, (r ^ a+), ^ ' 

and let us choose the convenient initial conditions i^(r,t)|j^o = 0, a{Vl,t)\t=Q = 0, (i7 = {O^iji) — the spherical angles), 
moreover (1 -I- /i(r))|r=Q = f2{f')\r=a (continuity condition), F{r,t)\r=a = 0^ / p{r,t)d?'r = (neutrality condition). 
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We thus arrive at the exphcit form of the solutions of Eqs ([20l) and pT|) (as it is described in the Appendix 

Mr) = _fcx|±ie-fcT(--) for r < a, 



/2(r) = [kra - ^ (l - e'^'^-")] ^^^gf^, /or r > a, 



(23) 



where kx — .1 '22Uis£l — _e_ por time-dependent parts of electron fluctuations we find: 



/^l rn — ~-l n—1 



and 



oo Z 

a{n,t)^Y: E i^/m(f^) [■^si«(wo/t)e-*/"«(l-'5ii) + QimW'5iz] 

/ — I m— — l 

+ E E E ^'"'» ;^2_f2;+iw^ Yim(r?)we / dri^ji{kniri)sin(u'^it)e-*/'''>, 

l = lm=-ln=l p ^ ^ ' ,>! Q 



(25) 



where ji(^) = y j^Ii+i/2i£,) is the spherical Bessel function, Yimi^) is the spherical function, — ojy 
are the frequencies of electron volume free self-oscillations (volume plasmon frequencies), Xni are nodes of the Bessel 
function luqi — tOp^ ^'^^ frequencies of electron surface free self-oscillations (surface plasmon frequen- 

cies), and kni = Xni/a; to' — y'w^ — are the shifted frequencies for all modes due to damping. The coefficients Bim 

and Aimn are determined by the initial conditions. As we have assumed that 5p{r,t = 0) = 0, we get Bim — and 
Mmn = 0, except for Z = 1 in the former case (of Bim) which corresponds to homogeneous electric field excitation. 
This is described by the function Qim{t) in the general solution ([25l) . The function Qim{t) satisfies the equation: 

— m — + 7^ — m — + ^iQim{t) 

[E,{R,t)6„,^ + V2iE,iR,t)S„a + Ey{R,t)S , 



where Wi = luqi = (it is a dipole-type surface plasmon Mie frequency!^). Only this function contributes the 
dynamical response to the homogeneous electric field (for the assumed initial conditions). From the above it follows 
thus that local electron density (within RPA attitude) has the form: 

p{r,t)=po{r)+piir,t), (27) 

where the RPA equilibrium electron distribution (correcting the uniform distribution rie): 

Po(r) = ("H\+/f^^' "^"^ (28) 
[ nehv), /o'^ r > a, r a+ ^ ' 

and the nonequilibrium, of surface plasmon oscillation type for the homogeneous forcing field: 

{0, for r < a, 
E QirnimMn) for r>a, r^a + . (29) 
m— — 1 

In general, F{r,t) (volume plasmons) and a{il,t) (surface plasmons) contribute to plasmon oscillations. However, in 
the case homogeneous perturbation, only the surface I = 1 mode is excited. 

For plasmon oscillations given by Eq. (|29l) one can calculate the corresponding dipole. 



/■ , 47r o 

D{R,t) = e / (frrp{r,t) = —eqiR,t)a^, (30) 

where Qii{R,t) = yJ^q.^{R,t), Qi^i{R,t) = ^J^qy{R,t), Qio{R,t) = ^J^q^{R,t) and q{R,t) satisfies the 
equation (cf. Eq. (f26| ). 



^ 2d 



PT) 

q{R,t) = —E{R,t). (31) 
m 
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III. LORENTZ FRICTION FOR NANOSPHERE PLASMONS 



Considering the nanosphere plasmons induced by the homogeneous electric field, as described in the above para- 
graph, one can note that these plasmons are themselves a source of the e-m radiation. This radiation takes away 
the energy of plasmons resulting in their damping, which can be described as the Lorentz friction^. This e-m wave 
emission causes electron friction which can be described as the additional electric field^. 



E, 



2 d^D{t) 



(32) 



where w = is the light velocity in the dielectric medium, and D{t) the dipole of the nanosphere. According to Eq. 
([30|) we arrive at the following: 



Substituting it into Eq. ((3T|) we get 



2_ a_ 

To dt 



q{R,t) 



'-E{R, t) 



3aji V V 



If rewrite the above equation (for E—Q) in the form 

q{R,t) = 



9^ 



2 

To 



2 
3a; 1 



thus the zeroth order approximation (neglecting attenuation) corresponds to the equation; 

yQ2 



q{R,t) = 0. 



(33) 



(34) 



(35) 



(36) 



In order to solve Eq. (j35p in the next step of perturbation, in the right-hand-side of this equation one can subsitute 

by -iofq{t) (acc. to Eq. 

Therefore, if one assumes the above estimation, 
renormalised damping term; 



9*3 



,2dq{t) 



then one can include the Lorentz friction in a 



2d_ 

Tdt 



q{R,t) 



where 



1 
r 



1 

tq 



2A 



B 



2a 



'■-E{R,t), 



(37) 



(38) 



(Ab is the free path in bulk, vp the Fermi velocity, and C ~ 1 a constant )^i^ 



where we used for — ~ t^^^ , „ 

To ^ 2Xb 2a _ _ 

which corresponds to inclusion of plasmon damping due to electron scattering on other electrons and on nanoparticle 
boundary. The renormalised damping causes the change in the shift of self-frequencies of free surface plasmons. 



Using Eq. psp one can determine the radius oq corresponding to minimal damping. 



ao = — (fFC^\/i/2) 



1/4 



(39) 



For nanoparticles of gold, silver and copper in air, in water and in a colloidal solution, one can find ao < lOnm 
(cf. Tab. 1), which corresponds to the experimental data^^— . For a > oq damping increases due to Lorentz 
friction (proportional to a^) but for a < damping due to electron scattering dominates and causes also damping 
enhancement (with lowering a, as ~ ^, cf. Fig. 1), 
which agrees with experimental observation a^'^'^^ . 
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FIG. 1: Effective damping ratio for surface plasmon oscillations, Eqs H37I) . (|38p . the upper (blue) curve is the sum of both 
terms, ~ ^ (red) and ~ (green); the minimum corresponds to minimal damping for radius ao, Eq. (|39p . left — for Ag in the 
air, right — for An in colloidal water solution 



Tab. 1. gq — nanosphere radius corresponding to minimal damping 



refraction rate of the surrounding medium, uq 


Au, flo [nm] 


Ag, ao [nm] 


Cu, ao [nm] 


(air) 1 


8.8 


8.44 


8.46 


(water) 1.4 


9.14 


9.18 


9.20 


(colloidal solution) 2 


9.99 


10.04 


10.04 



Surface plasmon oscillations cause attenuation of the incident e-m radiation where the maximum of attenuation is 
at the resonant frequency?^ cji = ^ojf — This frequency diminishes with rise of a, for a > ao according to Eq. 
(pS]). which agrees with experimental observations for Au and Ag presented in Fig. 2, and Tab. 2 (Au) and Tab 3 
(Ag). 



Tab. 2. Resonant frequency for e-m wave attenuation in Au nanospheres 



radius of nanosheres [nm] 


10 


15 


20 


25 


30 


40 


50 


hio'i (experiment) [eV] 


2.371 


2.362 


2.357 


2.340 


2.316 


2.248 


2.172 


huj[ (theory) [eV], no = 1.4 


3.721 


3.720 


3.716 


2.702 


3.666 


3.415 


2.374 


hu'i (theory) [eV], no = 2 


2.604 


2.603 


2.600 


2.590 


2.565 


2.388 


1.656 



Tab. 3. Resonant frequency for e-m wave attenuation in Ag nanospheres 



radius of nanosheres [nm] 


10 


20 


30 


40 


tioj'i (experiment) [eV] 


3.024 


2.911 


2.633 


2.385 


huj[ (theory) [eV], no = 1.4 


3.707 


3.702 


3.654 


3.410 


huj'i (theory) [eV], no = 2 


2.595 


2.591 


2.557 


2.384 



IV. PLASMON-MEDIATED ENERGY TRANSFER THROUGH A CHAIN OF METALLIC 

NANOSPHERES 

Let us consider a linear chain of metallic nanospheres with radii a in a dielectric medium with dielectric constant e. 
We assume that spheres are located along z-axis direction equidistantly with separation of sphere centers d > 2aM^. 
At time t = we assume the excitation of plasmon oscillation via a Dirac delta ~ 6{t) shape signal of electric field. 
Taking into account the mutual interaction of induced surface plasmons on the spheres via the radiation of dipole 
oscillations, we aim to determine the stationary state of the whole infinite chain. For separation d much shorter than 
the wavelength A of the e-m wave corresponding to surface plasmon self-frequency, the dipole type plasmon radiation 
can be treated within near-field regime, at least for nearest neighbouring spheres. In the near-field region a < Rq < X, 
the radiation of the dipole D{t) is not a planar wave (as for far-field region, Rq ^ A) but of only electric field type 
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FIG. 2: Extinction spectra for nanospheres of Au (a) and Ag (b) in colloidal water solution for various sphere radii 



in retarded form (without magnetic field)^: 



E{R, i?o, t) 



3n{ n- D i RA 



Rq 

V 



D\RA 



Ro 

V 



(40) 



R — position of the sphere (center) irradiating e-m energy due to its dipole surface plasmon oscillations, Rq position 
of another sphere (center), with respect to the center of the former one, where the field E{R,Ro,t) is given by the 
above formula, Rq < X, n = Rq/Rq, v = c/^/e = c/uq. 

When both vectors R and Rq are along the z-axis (the linear chain) the above equation can be resolved as: 



Ea{R, -Ro, t) 



Da{R,t^ 



Ro 

V 



(41) 



where a = {x,y,z), = cry = — 1 and az — 2. Assuming that the z-axis origin coincides with the center of one 
sphere in the chain, for the sphere located in the point Ri = (0,0, /d), an electric field caused by neighbouring 
spheres, E{Rm, Rmi,t), and the Lorentz friction force caused by self-radiation, EL{Ri,t), has to be considered. By 
virtue of Eq. (|3ip the equation for surface plasmon oscillation of the sphere is 



5t2 



2^d_ 
To dt 



q{Rl^t)^^ J2 
m ^ — ' 



E{Rm, Rml,t) H -EL{Ri,t), 



m— — oo, m^l^ Rmi<\ 



m 



(42) 
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provided that the dipole field of the m*^ sphere can be treated as homogeneous over the Z*^ sphere and the sum over 
m is confined by the distance of m*^ sphere from l^^ sphere not exceeding the near-field range (~ A). In the case of 
the equidistant chain, Ri = Id and Rmi — \l ^ m\d, and using Eqs (|4T|) . (pO| and [33]) . one can rewrite Eq. (|42l) in the 
form: 



dt'^ To dt 3uj 



2 /a;ia\3 



qa{ld,t) = a a 



m—oo 



d3 ^ 

m— — oo, m^/, \l—m\d<X 



Qa {md,t~ m\) 



(43) 



here a = x,y, which describe the transversal plasmon modes and a = z, which describes the longitudinal one. The 

' " ' 47r 
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above equation coincides with the appropriate one from Refs |23ll25j . if one assumes that ^a^rie = N = 1 and neglects 



the retardation of the field. 

Taking into account the periodicity of the infinite chain, one can consider the solution of the above equation in the 
form 



q{ld,t)=q{k,t)e-'''''^. 
The right-hand-side term in Eq. (|43l) attains the form 



(44) 



m— — oo, m^l 



qc{rnd.t~^\l-7n\) q^(md,t~ ^\l-m\) 



|/-m|3 



m— — oo 



|;-m|3 



Thus the Eq. (^5)) can be written as follows: 



1^ |/-m|3 
m=l + l 



2d 2 /wia\3 



To 9i 3a;i V w / 9t3 



qa{k,t) ^ Gal^l-^ ^ 



cos{mkd) ^ 

^— — -q{k,t- md/v) (45) 



m— 1, md<X 



This equation is linear and therefore we look for the solutions of the shape: qa{k,t) = Qa{k)e^'^°'* , and 

2iuJa 



■^1 + 



I - 2 



(46) 



where 



(i3 — ' rn 

m— 1, 7nd<X 



cosimkd) ( ujfyjnd 
T, cos ' 



(47) 



and 



1 



1 uj^a I u)na,\ 2 

TaiuJa) To ' 3 



— [ 1 + CTaWi 



E 



COS 



(mkd) sin (^^) 



rf3 — ' 

m— 1, md<X 



(48) 



If we confine the sum in the Eq. (147^ to m = 1 (the nearest neighbour approximation) we get 



cos{kd)cos 

6(3 \ i; 



(49) 



and from Eq. (|48 



2„3 



1 ujfa' 



T„(Wa) To 4u(i2 



— j ~ikd-n) +- 



, /or a = x,y 



(50) 



and 



1 



Tz(a;z) To 2?;(i2 



+ {kd-Tif 



, /or q; = z. 



(51) 
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In the derivation of two above formulae the following summation was performed^S; 

_2_ j2 cos(mfcd) ^^^ f^uj,^^ ^ _2_ J2 :^[sin{km,d + ujamd/v) ~ sin{kmd- ujamd/v)] 
--^kd 



4 ^tP~ 



as the terms in the sum drop quickly to zero then the above formula well approximates the sum with limitation 
md < A. 

Assuming now cj^ = + iuj^ the Eq. (|46p gives the dependence of uj^ and cj^ on fc. The general solution of Eq. 
(|43| attains the form, 



qa{ld,t) = ^ Qa{kn 



J(u: ^(k„)t-knld)-LO ^ {k„)t 



(52) 



where fc„ — L — Ngd is assumed length of the chain with Ng spheres, and periodic (of Born-Karman type) 

boundary condition imposed. The components of Eq. (|5^ describe monochromatic waves with wavelength A„ — 
1^ = ^, which are analogous to planar waves in crystals, when damping is not big, i.e., when uj^ <^ uj^. Provided 
this inequality one can approximate: 
for transversal modes (a = x, y) 



1 1 ula? 

Ta To 4ua'^ 

and for longitudinal mode (a = z) 

"11 w?a3 

UJ= — = 



2a^ 

1 H — —cos(kd)cos(uj d/v) 
d-^ 



(53) 
(54) 



1 - 



4a^ 
^3" 



■cos{kd)cos{uj^d/v) 



Tz To 2vd'^ 



+ {kd-T:Y - 



(55) 
(56) 



From the Eqs ([5l| and ([5S1) it follows that can change its sign. In the case of < the oscillations are 
destabilized, which could be avoided by inclusions of some nonlinear terms neglected in the expression for the Lorentz 
friction, which in more accurate form^-3. includes also a small nonlinear term with respect to D, aside from the term with 
^Q§-- Including of it will result in damping of too highly rising oscillations leading to stable amplitude of oscillations. 
Due to this stabilisation caused by nonlinear effects, undamped wave modes of dipole oscillations will propagate in 
the chain in the region of parameters where uj^ < Q (and with fixed amplitude accommodated by nonlinear term). 
The condition aj^, = = 0, for critical parameters, resolves into: 



iiJad 



{kd 



TT 



for a — {x, y) and for a — z, 



^\ =-{kd-ny + — 



Avd^ 
TQUjfa-^ 



2vd^ 



(57) 



(58) 



Obtained from the above equations (^^^) leads to determination of the dependence of wave vector k with respect to 
parameters d and a, via Eqs (|53p -(|56 p . Solution for this equations, found numerically for the chain of Ag nanospheres, 
is depicted in the Fig. 3. 

The undamped plasmon waves in the chain appear if d < dmax and have k — Tr/d, dmax — 98.5(68.8) nm for 
transversal(longitudinal) modes. For example, for Ag spheres with radius a = 20nm and separation d = 60 nm, the 
undamped transversal modes appear for < fed < 7r/4 or 37r/4 < kd < 2tt and longitudinal for 37r/4 < kd < 57r/4. 

Let us underline that the determined undamped plasmon oscillation wave modes explain the numerically observed 
similar behaviour—. 
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FIG. 3: Wave vector kd versus sphere separation d, at constant d/a = 3 (a — sphere radius) (upper), for a = 20 nm (lower) for 
Ag sphere chain, transversal modes — red, longitudinal mode — blue 



V. CONCLUSIONS 



In the present paper we analysed plasmons in large metallic nanospheres induced by homogeneous time-dependent 
electric field. Within all- analytical RPA quasiclassical approach the volume and surface plasmons are described and a 
proof that only dipole-type of surface plasmons can be induced by a homogeneous field (while none of volume modes) 
is given. An irradiation of energy by plasmon oscillations is described within the Lorcntz friction effect. Its scaling 
with the nanosphere dimension leads to sphere radius dependent shift of resonant frequency, similarly as observed in 
experiments. The description of surface dipole-type plasmon oscillations in single nanospheres is applied to analysis of 
collective oscillation in linear chain of metallic nanospheres. The wave-type collective plasmon oscillations in the chain 
are also considered. The undamped region of wave energy transport through the chain is found for a certain sphere 
separation in the chain with corresponding appropriate wavelength of plasmon waves. This phenomenon confirms a 
similar behaviour observed by numerical simulations^. 
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Appendix A: Analytical solution of plasmon equations for the nanosphere 



Let us solve first the Eq. ([^0]) . assuming a solution in the form: 

Spi {r,t) = He [fi (r) + F{r, t)] , for r < a, 



(Al) 
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Eq. (|20|) resolves thus into: 

'■F(r,t) 2_dF{r 
0*2 ^ TO at 

The sohition for function (nonsingular at r = 0) has thus the form 



(A2) 



(A3) 



where a — const., kx = \f^^f^ ~ ^^'^ — inverse Thomas- Fermi radius), ujp — y^^^^ (bulk plasmon 

frequency) . 

Since we assume F{r,0) = 0, then for function F{r,t) the solution can be taken as, 

F{r,t)^F^ (r)sm(a;'t)e"^°* (A4) 



where w' = a/w^ + I/tq . F^{r) satisfies the equation (Helmholtz equation): 

V^FUr) + eF^{r)=0, (A5) 

with A:^ — — TTT- ^ solution of the above equation, nonsingular at r = 0, is as follows: 

FUr) ^ Aji{kr)Yiram, (A6) 

where A — constant, — y/''^ / {'^£,)h+i/2i.i) — the spherical Bessel function [In{£,) — the Bessel function of the first 
order], Yim{^) — the spherical function (f2 — the spherical angle). Owing to the boundary condition, F(r,t)|r=a = 0, 
one has to demand ji{ka) = 0, which leads to the discrete values of fc = k^i = Xni/a, (where Xni, n- — 1,2,3..., are 
nodes of ji), and next to the discretisation of self- frequencies oj: 

^li = (l + ^) • (A7) 
The general solution for F{r^t) has thus the form 

oo I oo 

F{r,t)^Yl $^^/W/(fcnir)yi™(l])sin(w;it)e-"°*. (A8) 

l—Q m——l n—1 

A solution of Eq. (f2T|) we represent as: 

Sp2{r, t) = ne/2(r) + a{fl, t)6{r + e-a), for r > a, {r a+, i.e.e 0). (A9) 
The neutrality condition, / p{r, t)d^r = N^, with Sp2{r, t) = cr(a;, t)6{a + e — r) + rie f2{r) , { e — 0), can be rewritten 

a oo a 

as follows: ~ j drr^ fi{r) = J drr^ f2{r), J d'^rF{r,t) — 0, J dil(T{il.,t) = 0. Taking into account also the continuity 

a 

condition on the spherical particle surface, 1 -t- /i(a) = /2(a): one can obtain: f2{r) — /3e~'^^^''~°-'/(^T'') and it is 
possible to fit a (cf. Eq. dM])) and (3 constants: a = j3 = kTa- -^^^ifii (l - e"2feTa)_which gives Eqs 

a 

From the condition J d^rF{r,t) = and from Eq. (|A8|) it follows that ^oon — 0, (because of / d^Yim{(jj) = 


47r(5/o^mo)- 

CO 

To remove the Dirac delta functions we integrate both sides of the Eq. (l^lT) with respect to the radius (J r^dr...) 



and then we take the limit to the sphere surface, e — > 0. It results in the following equation for surface plasmons: 

^^ + ^^ = -£ E u:ln^{n)jdnM^,^t)Y*^{n,) 

1=0 m,=-l 

00 I ao a ^ij_2 fAlO) 

-t-w^ne X] J2 ^inin2r+i^hn{^) J dri-^ji{kniri)sin{ujnit), 

+ ^^4^3 [E,{R, t)Yiom + V2E,iR, t)Yn i^) + V2Ey{R, t)Yi^i{n)] , 
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where Wqj = Wp^iqrr ■ In. derivation of the above equation the following formulae were exploited, (for a < ri): 



d 1 d ^ a', ^ ^ /a'-i 



da + rf - larxcos^ da ^ ^ r[ 



J2 —Piicosj) = J2 -^Pii^os^). (All) 



where Pi{cos^) is the Legendre polynomial [Pi{cos^) — ^ Yimi^)Yi^{Qi)], 7 is an angle between vectors 

m— — l 

a = ar and ri, and (for a > ri): 

Taking into account the spherical symmetry, one can assume the solution of the Eq. (jAlOp in the form: 

00 / 

a(17,t) = ^ ^ g,„,(i)y,„(r!). (A13) 

1=0 m=-l 

From the condition J a{ujt)dfl = it follows that goo = 0. Taking into account the initial condition cr(w, 0) = we 
get (for I > 1), 

qim{t) = ^sin{Lu'^,t)e-'/^"{l - Sn) + Qi„,{t)Sn 



where = \/^'oi ~ -'-/''"o ^^^^ Qim{t) satisfies the equation: 
d'^Qi,n{t , 2 dQim{t) , 2 ^ u\ ^ 



Thus cr(uj,t) attains the form: 



+ ujl,Qi.m{t)^ —^/^ E,{R,t)6^o + V2E,{R,t)6„^i+V2Ey{R,t)6^_i . (A15) 



00 I 1 

a(l),t) = E E i^i™(f^)^sin(w^,,Oe"*/"" + E Qim(imm(f^) 

I— 2m— — I rn——l 



(AI6) 

'"i ;w2-(2/+i)l. , 



+ E E E ^»i"t ;^2_(2;+iw^ Y;„(r?)ne/dri^j;(fc«iri)sm(^;;Oe~*^^°- 



^ S. Pillai, K. B. Catchpole, T. Trupke, G. Zhang, J, Zhao, and M. A. Green, Appl. Phys. Let., 88, 161102 (2006) 

^ M. Westphalen, U. Kreibig, J. Rostalski, H. Liith, and D. Meissner, Sol. Energy Mater. Sol. Cells 61, 97 (2000), M. Gratzel, 

J. Photochem. Photobiol. C: Photochem. Rev. 4, 145 (2003) 
^ H. R. Stuart and D. G. Hall, Appl. Phys. Lett. 73, 3815 (1998); H. R. Stuart and D. G. Hall, Phys. Rev. Lett. 80, 5663 

(1998); H. R. Stuart and D. G. Hall, Appl. Phys. Lett. 69, 2327 (1996) 

* D. M. Schaadt, B. Feng, and E. T. Yu, Appl. Phys. Lett. 86, 063106 (2005) 

5 K. Okamoto, et al.. Nature Mat. 3, 661 (2004); K. Okamoto, et al., Appl. Phys. Lett. 87, 071102 (2005) 
" C. Wen, K. Ishikawa, M. Kishima, K. Yamada, Sol. Cells 61, 339 (2000) 
L. Lalanne, J. P. Hugonin, Nature Phys. 2, 551 (2006) 

* A.V. Zayats, I. I. Smolyaninov, and A. A. Maradudin, Phys. Rep. 408, 131 (2005) 
® S.A. Mayer, Plasmonics: Fundamentals and Applications, Springer VL 2007 

^° W. L. Barnes, A. Dereux, and T. W. Ebbesen, Nature 424, 824 (2003) 
" N. Engeta, A. Salandriw, and A. Alu, Phys. Rev. Lett., 95, 095504 (2005) 
^2 S. A. Maier and H. A. Atwater, J. Appl. Phys., 98, 011101 (2005) 

G. Mie, Ann. Phys. 25, 329 (1908) 
" M. Brack, Phys. Rev. B 39, 3533 (1989) 
1^ L. Serra et al. , Phys. Rev. B 41, 3434 (1990) 

M. Brack, Rev. of Mod. Phys. 65, 677 (1993); 

W. Ekardt, Phys. Rev. Lett. 52, 1925 (1984) 



14 



C. F. Bohrcn, D.R. Huffman, Absorption and Scattering of Light by Small Particles, Wiley, New York (1983); U. Kreibig, 
M. VoUmer, Optical Properties of Metal Clusters, Springer, Berlin (1995); J. I. Petrov, Physics of Small Particles, Nauka, 
Moscow (1984); C. Burda, X. Chen, R. Narayanan, M. El-Sayed, Chem. Rev. 105, 1025 (2005) 

^® D. Pines, Elementary Excitations in Solids, ABP Perseus Books, Massachusetts (1999) 
^° L. Jacak, J. Krasnyj, A. Chepok, Fizika Niskich Temp. 33, (2009) 

D. Pines and D. Bohm, Phys. Rev. 85, 338 (1952); D. Bohm and D. Pines, Phys. Rev. 92, 609 (1953) 
L. D. Landau and E. M. Lifshitz, Field Theory, Nauka, Moscow (1973) (in Russian) 

M. L. Brongcrsma, J. W. Hartman, and H. A. Atwater, Phys. Rev. B 62, R16356 (2000) 
2* U. Kriebig and L. Geinzel, Surf. Sci., 156, 678, (1985) 

2"'' S. A. Maier, P. G. Kik, and H. A. Atwater, Phys. Rev. B, 67, 205402 (2003) 

F. Stietz, I. Bosbach, T. Wenzel, T. Vartanyan, and A. Goldmann, F. Trager, Phys. Rev. Lett., 84, 5644 (2000) 

F. Stietz et al, Phys. Rev. Lett. 84, 5644 (2000) 
^® M. Scharte et all., Appl. Phys. B: Laser Opt. 73, 305 (2001) 
^® I. S. Gradstein, I. M. Rizik, Tables of Integrals, Fizmatizdat, Moscow (1962). 
^° V. A. Markel and A. K. Sarychcv, Phys. Rev. B, 75, 085426 (2007) 
3^ E. Hao, R. C. Bayley, G. C. Schatz, J. T. Hupp, and S. Li, Nano Lett. 4, 327 (2004) 

B. Lamprecht, A. Leitner, and F. R. Aussenegg, Appl. Phys. B: Lasers Opt. 64, 269 (1997) 



